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ABSTRACT 

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster 
Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and 
SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED). We present maps of dust temperature, dust 
mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature 
in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is 
symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with 
a large sample of resolved galaxies to be observed by Herschel. 
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1. Introduction 

Infrared data have been widely used to determine the composi- 
tion and distribution of dust in galaxies since the launch of the 
Infrared Astronomical Satellite (IRAS) in the 1980s. However, 
dust masses for nearby galaxies calculated from IRAS 60 and 
100 jim measurements were found to be a factor of ten lower 
than expected w hen compared to the M ilky Way gas-to-dust ra- 
tio of 100-200 ( Devereux & Young|[T99 0). The Milky Way dust 
mass was calculated by measuring the depletion of metals from 
the gaseous phase of the interstellar medium (ISM) and by com- 
paring gas column densities to dust extinction (Whittet 20031 
its value implying that most of the dust- mass emits radiation 
at wa velengths longer than 100 jim (e.g., Devereu x & Youn g 
Il990h . Analyses of data from the Spitzer Space Tele scope hav e 
determined gas-to-dust ratios of -150 (e.g., Draine et al.ll2007l) . 
However, these analyses were generally limited by the number, 
and low signal-to-noise ratio of data at wavelengths greater than 
160 jim. These studies therefore had difficulty detecting emis- 
sion from dust with temperatures lower than 15 K, and their re- 
sults were biased towards warmer dust temperatures and lower 
masses. Owing to the difficulties at these wavelengths, high res- 
olution studies of galaxies have been previously limited. These 
studies are important for understanding how dust interacts with 
the other phases of the ISM, the sources of dust heating and how 
the distribution and temperature of dust varies with mo rphology. 
The Herschel Space Telescope (JPilbratt et all 12010 ), with its 
two photometric instruments, constrains both sides of the peak 
in the far infrared spectral energy distribution (SED, see Fig. 
[T]). The Photodetector Array Camera & Spectrometer (PACS, 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



iPoglitsch et all 1201 Oh has three photometric bands at 70, 100, 
and 160 jim at superior angular resolutions to those provided 
by Spitze r. The Spectral and Photometric Imaging Receiver 
(SPIRE, Griffin et al 2010) also has three photometric bands ob- 
serving simultaneously at 250, 350, and 500 jim with high sen- 
sitivity and angular resolution. 

The Herschel Virgo Cluster Survey (HeViCSQ) is a Herschel 
Open Time Key Project that will observe -64 deg^ of the Virgo 
cluster. This will provide a large sample of resolved galaxies be- 
cause about 48 late-type galaxies will be observed by HeViCS 
with optical diameters larger than 3'. In this Letter, we present 
an insight into what will be possible with the full HeViCS sur- 
vey by applying a resolved dust analysis to infer dust temper- 
atures, surface densities, and gas-to-dust ratios for NGC 4501 
and NGC 4567/8. These galaxies were chosen because they are 
among the largest angular size systems in the HeViCS Science 
Demonstration Phase (SDP) field. In Sect. 2, we present the ob- 
servations and data reduction, and in Sects. 3 and 4, we present 
our analysis and results, respectively. 

2. Observations and data reduction 

As part of the Herschel SDP, a 4x4 sq deg field of the Virgo 
cluster centred approximately on M87 has been observed in par- 
allel mode, simultaneously at 100 and 160 jim with PACS and 
at 250, 350, and 500 jim with SPIRE. For details of the obsery - 
ing strategy and data reduction, we refer to lDavies et all (1201 Oh . 
Following the recommendations given by the PACS and SPIRE 
instrument control centres (ICCs), the 100 and 160 jim PACS 
flux densities were scaled by dividing by 1.06 and 1.29, respec- 
tively, and SPIRE flux densities by multiplying by 1.02, 1.05, 



More details on HeViCS can be found at |http ://w w w.hevic s . org | 
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Fig. 1 : Global and example single pixel one component SED fits. 
If not shown, uncertainties are smaller than the symbol size. The 
70 jim point (blue) is used as an upper limit (see SectO. 



3. Analysis 

We fit simple models to the SEDs of each pixel in our galax- 
ies to produce maps of estimated dust temperature and dust 
mass. All images were first convolved to the resolution of the 
500 jim image (which has a PSF of the largest FWHM), us- 
ing a customised ke rnel created using the procedure outlined in 
iBendo et aP (l2010bh . The images were then regridded to the 14'' 
pixel size of the 500 jim map. We note that since the pixel size is 
smaller than the 500 jim beam size, neighbouring pixels are not 
independent. Owing to lower sensitivities in the PACS 100 jim 
band and negative artifacts created around the galaxy during 
the current PACS data reduction process, we only consider pix- 
els with a flux density >10cr in the convolved and regridded 
100 jim image; the estimated sensitivities are 7.3 mJy pix"^ for 
the cross-scan region and 10.4 mJy pix"^ for the single scan 
(see lDavies et alll201Ql for optimal resolution sensitivities). Thus, 
in the present work we do not study the outermost regions of 
the galaxies; the full galaxy extent will be considered in future 
works. 

We fit the SEDs in the 70-500 jim range with one and two 
component modified black-body models. The equation for a one 
component black body is given by 



and 0.94 for the 250, 350, and 500 jim bands, respectively. The 
point spread function (PSF) of the PACS images has a FWHM of 
12.70'' X 6.98" and 15.65" x 11.64" in the 100 jim and 160 jim 
band, respectively (where the orientation depends on scan direc- 
tion). The SPIRE PSF has a FWHM of 18.1", 25.2", and 36.9" 
at 250 jim, 350 jim, and 500 jim, respectively. The calibration er- 
ror is assumed to be 20% for the PACS band and 15% for the 
SPIRE bands. The PACS and SPIRE images are shown in Figs. 
[2]and[3l 

This Letter focuses on NGC 4501, which is an Sbc galaxy 
and has an optical diameter (major axis) of 7.23', and the galaxy 
pair NGC 4567/8, which are both of Sc type with optical diam- 
eters of 2.92' and 5.1', re spectively (values taken from GOLD 
Mine: lGavazzi et al.ll2QQ3i) . As NGC 4501 is on the edge of the 
SDP field, observations were only obtained in one scan direc- 
tion. 

In addi tion to the Her s chel d a ta, we used S p itzer 70 jim 
data from iKennicutt et"aD (l2003h . iKennev et all (l2006b . and 
IStrucket aP (l2004l), which were repro cessed using the tech- 
niques described in lBendo et al.l (1201 Obi) . The calibration uncer - 
tainty were assumed to be 5% at 70 jim (Gord on et al.ll2007l) . 
We also use GALEX FUV {A = 1539A,A^ = 442A) data 
downloaded from the MAST (G4/5) archive. For both NGC 
4501 and NGC 4567/8 , we used Hi maps from the VIVA 
Hi survey (IChung et al.1 12009), and for NGC 4501 we also 
used CQ(/= 1-0) maps from the Nobeyama CO ATLAS sur- 
vey (iKunoetal. 1 120071 The Hi and CO images are shown in 
Figs. |2] and [31 For NGC 4567/8, w e used global CO(/=1-0) 
measurements from lono et al. ('2005'). In the following analysis, 
neutral atomic hydrogen masses were calculated using Mm = 
2.36 X lO^D^p^F/ (Mo), where FI is the integrated HI line flux 

in units of Jy kms"^ ( lRobertsiri962t) and Dmpc is the distance in 
Mpc. We assumed the distan ce to the c entre of the M87 cloud 
in the Virgo cluster given bv lMeietafl (120071) of 16.7 Mpc. To 
convert a CO integrated line intensity to a value of molecular 
hydrogen column density, we as sume an X factor of 1.9 x 10^^ 
cm"^[K km s"^]"^ derived from 'Strong & Mattox| (Il996l) using 
models of gamma ray scattering. However, we note that the value 
of the X factor is notoriously uncertain an d may also vary, for 
example with metallicity (^e.g. Jlsraelll2000h . 



Fy = ^MBAT). 



(1) 



Here M is the dust mass, T is the dust temperature, By(T) is the 
Plank function, D is the distance to the galaxy (see Sect.[2l), and 
Ky is the dust emissivity. The dust emissivity is assumed to be a 
power law in this spectral range, where Ky oc v^. We assume a 
value of /Co = 0.192 m^kg-^ at 350 jim (Draine 2003) and/5 = 2, 
though we note that the values of k andyS are notoriously uncer- 
tain. The best-fit solution is found by minimising the chi- squared 
(X^) function. To adjust for the filter-band passes, the SED is 
convolved with the filter transmission in the fitting process. 

We find that the data in the 100-500 jim range is accurately 
fitted by a single black-body component. However, for the ma- 
jority of pixels a significant flux excess is measured at 70 jim 
relative to this fit. This indicates that a warmer dust compo- 
nent is pres ent, which is consistent with the results found by 
iBendo et aP ([2010a) . A two-component model provides a good 
fit by using the same component through the 100-500 jim data 
as the one-component model and adding a warm component that 
fits the 70 jim emission. The warm component, however, is not 
tightly constrained since a wide range of warm temperatures sta- 
tistically provide equally good fits, due mainly to the small num- 
ber of data points on the Wien side of the SED. Therefore, in this 
work we choose to use a one-component model with the 70 jim 
data point as an upper limit. This choice has little eff'ect on the to- 
tal dust mass since the cold component dominates the dust mass. 
Example global and single-pixel SEDs are shown in Fig. [TJ We 
note that for most of the SED fits a ~ Icr excess is found for the 
250 jim band, for which we do not find any immediate explana- 
tion and we plan future investigations. 

In addition to maps of dust temperature, we also produced 
maps of dust mass surface density and the gas-to-dust ratio. For 
NGC 4501, we calculated Mgas = Mm + Mh^; for NGC 4567/8, 
we do not have a CO(/=1-0) map, so we only consider the Hi 
component of the gas mass. 

The uncertainties in the flux densities are calculated by 
adding in quadrature the calibration uncertainty, background 
subtraction uncertainty, and the uncertainty from the error map 
(created by the standard pipeline). For the pixels used in this 
analysis, the calibration uncertainty dominates for the PACS and 
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Fig. 2: Images of NGC 4501 in their native resolution and pixel size. The first three columns show the MIPS 70 jim PACS, SPIRE, 
and gas images on a log colour scale. The right-hand column shows maps of dust temperature, dust surface density (log scale), and 
gas-to-dust ratio. Beam sizes for the respective wavelengths are shown. The panels are 6'x6' in size and centred on 12h32m00.0s, 
-hl4d25ml2.0s. 



SPIRE bands, except for the 500 jim band where the total un- 
certainty ranges from 15-23%. The average ;^^ value for all the 
fits is 2.5 with a standard deviation of 1.0, where the 90 % confi- 
dence limit is 6.3 for 3 degrees of freedom. To obtain an estimate 
of the uncertainty in the fitted parameters, a bootstrap technique 
is used. This yields an average uncertainty in the pixel-by-pixel 
fitting of 1.0 K for temperature and 23 % in the dust surface 
density. 



4. Results 

Maps of the cold dust temperature, dust-mass surface density, 
and gas-to-dust ratio are shown in Figs. [2] and [3l for NGC 4501 
and NGC 4567/8, respectively. The extent of the maps corre- 
sponds to ~ 65 % of the optical radius. In each case, we find 
a centrally peaked, symmetrical distribution of dust mass in the 
galaxies. However, the dust temperature does not have the same 
symmetric distribution and varies across both galaxies, peak- 
ing at ~ 22 K and decreasing to ~ 19 K towards the outskirts 
of the galaxies. For NGC 4501, the peak of the dust tempera- 
ture distribution corresponds to the bright spiral arm region seen 
to the southeast of the galaxy in the 70-160 jim images. For 
NGC 4567/8, the dust temperature distribution is also asymmet- 
ric, with higher temperatures in NGC 4568 that are coincident 
with peaks in the GALEX FUV. For both NGC 4501 and NGC 
4567/8, the global temperature is 20 K (and also for NGC 4567 
and NGC 4568, individually; we included an additional uncer- 
tainty in the measurements to account for the overlap in emission 



from the two galaxies). In the regions that we analysed, we saw 
no evidence of any excess of flux at 500 jim that would imply an 
even colder dust component is present. 

The dust-mass surface density distribution for both NGC 
4501 and NGC 4568 has a smooth gradient, decreasing by 
a factor ~ 4 from higher values at the centre to lower val- 
ues in the outskirts. This i s consistent with results found by 
iMunoz-Mateos et aP (l2009l) . who find exponentially decreas- 
ing dust surface densities for the SINGS galaxy sample. For 
NGC 4501 and NGC 4567/8, we find global dust masses of 
(1.2 ± 0.3) X 10^ Mo and (8.6 ± 1.8) x 10^ M©, respectively 
((1.7 ± 0.4) X 10^ Mo and (6.6 ± 1.5) x 10^ Mq for NGC 4567 
and NGC 4568 individually, respectively). 

For NGC 4501, the distribution of the gas-to-dust ratio is 
asymmetric, peaking in the southeast and varying by a factor ~2 
across the galaxy. The average gas-to-dust ratio for the pixel- 
by-pixel analysis is 52+13 for NGC 4501, which is in good 
agreement with the global value of 51 + 14. The variation in the 
gas-to-dust ratio within NGC 4501 is consistent with gas-to-dust 
profiles within our fitted radius for other Sb galaxies presented 
in Munoz-Mateos et al. (2009). For NGC 4568, there is a large 
reduction in the Hi-to-dust ratio in the central regions, which 
is probably due to higher molecular gas densities. The global 
(neutral) gas-to-dust ratios for NGC 4567 and NGC 4568 are 
32+12 and 39+14, respectively. Though in all cases the global 
gas-to-dust ratios are low compared to the Galactic value, they 
are consistent with the lower end of the range of values in 
iDraine et aP (l2007b (correcting for the diff'erent X factor used). 
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Fig. 3: Same as Fig.[2]but for NGC 4567/8. The GALEX FUV image has been smoothed slightly for display purposes. The panels 
are 5.4'x5.4' in size and centred on 12h36m33.6s, +Ildl5m00.0s. 



In lMunoz-Mateos et al.l (l20Q9h . Sb and Sc galaxies have higher 
gas-to-dust ratios at larger radii than investigated in the current 
work. For NGC 4501 , the global value is in good agreement with 
the average pixel-by-pixel value, which suggests that a similar 
increase with radius may not be present for this galaxy. We note 
that s ince NGC 4501 and NGC 4568 are mildly Hi deficient (0.4 
- 0.6 jChung et al.ll2QQ9l) . low global gas-to-dust values of these 
Virgo galaxies could arise from the stripping of the out er gas and 
dust disk caused by environmental eff'ects ( Cortese e t al 2010). 
For the large sample of galaxies to be observed by HeViCS, we 
will be able investigate whether this is a general feature of the 
cluster environment. 



5. Conclusions 

We have fitted SEDs on a pixel-by-pixel basis to investigate the 
dust mass, temperature, and gas-to-dust ratios of NGC 4501 and 
NGC 4567/8. We have measured dust surface densities that peak 
at the centre of the galaxy and decrease towards the outer re- 
gions. In contrast, the temperature distribution is asymmetric, 
higher temperatures peaking at some distance away from the 
centre of the galaxy and then decreasing with increasing radius 
towards the outskirts. Once the full HeViCS observations are 
complete, we will be able to extend the analysis to larger radii 
and lower dust surface densities. For ~48 late-type galaxies with 
optical radii greater than 3' to be observed with HeViCS, we will 
be able to begin to tackle questions about the dominant source 
of dust heating, how the dust interacts with the other phases of 
the ISM, and how morphology influences the dust distribution 
and temperature. By combining the results from HeViCS and 



the Herschel Reference Survey (iBoselli et al .1120101) . we will be 
able to study how environment affects the dust distribution and 
other properties of a galaxy. 

Acknowledgements. Thanks to Edward Gomez for assistance with APLpy im- 
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